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Abstract 

In this paper, we investigate the energy problem in general rela- 
tivity using approximate Lie symmetry methods for differential equa- 
tions. This procedure is applied to Bardeen model (the regular black 
hole solution). Here we are forced to evaluate the third-order approx- 
imate symmetries of the orbital and geodesic equations. It is shown 
that energy must be re-scaled by some factor in the third-order ap- 
proximation. We discuss the insights of this re-scaling factor. 
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1 Introduction 

Energy-momentum is an important conserved quantity in any physical the- 
ory whose definition has been under investigation for a long time from the 
general relativity viewpoint. Unfortunately, there is still no generally ac- 
cepted definition of energy-momentum for gravitational fields. The main 
problem is with the expression defining the gravitational field energy part. 
To obtain a meaningful expression for energy, momentum and angular mo- 
mentum for a general relativistic system, Einstein proposed an expression. 
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called Einstein's energy-momentum complex [T] . After that many complexes 
have been found, for instance, Landau-Lifshitz [2], Tolman [3], Papapetrou 
f5], M0ller [5], Weinberg and Bergmann [7]. Some of these definitions are 
coordinate dependent while others are not. Also, some of these expressions 
cannot be used to define angular momentum. 

Virbhadra [S] gave marvellous idea about coincidence of some complexes 
which attracted many researchers. Aguireregabiria et al. [H] proved that 
Einstein, Landau-Lifshitz, Papapetrou, Weinberg prescriptions yielded the 
same distribution of energy for any metric of the Kerr-Schild class if calcu- 
lations are performed in Kerr-Schild coordinates. Vagenas [10] and Sharif 
[TT] proved that these four prescriptions coincide for the (2 + l)-dimensional 
rotating BTZ black hole (named after Banados, Teitelboim and Zanelli) and 
regular black hole solution respectively. Xulu [12] extended this investigation 
and found the same energy distribution for the dyadosphere of the Reissner- 
Nordstrom (RN) black hole. However, all these prescriptions have their own 
drawbacks. 

The concept of approximate symmetry theories comes from the combi- 
nation of Lie group theory and perturbations. Using this idea two different 
so-called "approximate symmetry" theories have been constructed. The first 
theory was proposed by Baikov, Gazizov and Ibragimov [13] and the second 
was presented by Fushchich and Shtelen [H]. It may be mentioned here 
that a manifold having no exact symmetry can possess approximate sym- 
metries. As compared to exact symmetry, one can acquire more interesting 
information from a slightly broken (approximate) symmetry. This idea of 
symmetries would be more accurate as it does not depend on pseudo-tensors 
and hence agrees with the principles of general relativity. In order to discuss 
the approach of approximate symmetries, we use the concept of Lie sym- 
metries which is substantiated in spacetime isometries as well as in Noether 
symmetries. 

Kara et al. [T5] used approximate Lie symmetry method [13] to discuss 
the conservation laws of the Schwarzschild spacetime. Later, using the same 
procedure, in the second-order approximation, the lost conservation laws 
for RN [TB], Kerr-Newmann [T7] and Kerr-Newmann AdS [IH] spacetimes 
were recovered and some energy re-scaling factors were obtained. In a recent 
paper [19], we have discussed re-scaling of energy in the stringy charged black 
hole solutions using approximate symmetries. In going from the Minkowski 
spacetime to non-flat regular black hole solution, the recovered conservation 
laws (recovered in first-order approximation) are lost. It is expected that by 
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minimizing the charge, one should recover all lost conservation laws. In this 
paper, we use the same procedure to formulate energy re-scaling factor for 
regular black hole solution. 

The format of the paper is as follows. In the next section, we review 
the basic material about approximate symmetry methods for the solution of 
differential equations (DEs) and also give exact symmetries of the Minkowski 
spacetime and the first-order approximate symmetries of the Schwarzschild 
spacetime. In section 3, we study approximate symmetries of the regular 
black hole solution. Finally, we summarize and discuss the results. 



2 Mathematical Formalism 

Symmetry is a point transformation (a transformation which maps one point 
{x,y) into another {x*,y*)) under which the form of DE does not change. 
These transformations form a group known as point transformation group. 
Symmetries are important for their direct connection with the conservation 
laws through the Noether theorem [20] • It has been mentioned that if, for 
a given system of DEs, there is a variational principle, then a continuous 
symmetry invariant under the action of the given functional provides corre- 
sponding conservation law j2T| [22]. 

A system of p ODEs, each of n^'^-order j23| [23] 

E„(5; x(s), x'(s), x"(s), x(")(s)) =0, (a = 1, 2, 3..,p) (1) 

(here x is dependent variable, s is independent variable and x', x",.--)X^"^ 
denote the first, second and so on ?7,*'*-order derivative of x w.r.t. s) under 
the point transformation: (s,x) — y {^{s,x),ri{s,x)) can admit a symmetry 
generator 

X = e(^,x)| + x)^ (2) 

if and only if, on the solution of the ODEs, = 0, the following symmetry 
condition 

xM(E„)|e..o=0 (3) 
is satisfied. Here n^'^-order extension of the symmetry generator given by 
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Eq.(|2]) can be written as 



+ ^:^„)(.,x,x',...,x("))^ (4) 



and the corresponding prolongation coefficients are 

The fc^'^-order approximate symmetries of a perturbed ODE 

E = Eo + eEi + e2E2 + ... + e'=Efc + 0(e'=+^); e E R+ (6) 
is given by the generator 

X = Xo + eXi + e2X2 + ... + e'=Xfc (7) 
if the following symmetry condition holds [25] 

= [(X = Xo + eXi + e2X2 + ... + e'=Xfc)(E = Eo + eEi 

+ e^E2 + ... + e^Efc)] |E=Eo+eEi+e2E2 + ...+e'=Efc= 0{t^^^). (8) 

Here Eq is the exact part of ODE, Ei, E2 are respectively the ffist and second 
order approximate parts of the perturbed ODE and so on. Xq represents 
the exact symmetry generator and Xi, X2 denote respectively the ffist and 
second order approximate parts of the symmetry generator and so on. For the 
/c*'^-order approximate symmetry generator, the terms involving e'-*^^^-' and its 
higher powers can be substituted to zero so that RHS of Eq.(IH]) becomes zero. 
As compared to exact symmetries, these approximate symmetries do not 
necessarily form a Lie algebra rather do form the so-called "approximate Lie 
algebra" (up to a specified degree of precision) ^26j. The perturbed equation 
always admits an approximate symmetry eXo called a trivial symmetry. If a 
symmetry generator X = Xq + eXi exists with Xq 7^ and Xi 7^ A;Xo, k is 
an arbitrary constant, then it is called non-trivial symmetry |27] . 

It has been mentioned [28] that 10 Killing vectors (Appendix A) of 
Minkowski spacetime provide conservation laws for energy and linear mo- 
mentum as well as for angular and spin angular momentum. The algebra 
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calculated from the geodesic equations contains some additional symmetries 
which provide no conservation law [22]. The Schwarzschild spacetime has four 
isometrics Yq, Yi, Y2, Y3 which correspond to conservation laws of energy 
and angular momentum only. The symmetry algebra of the Schwarzschild 
spacetime (calculated through the geodesic equations) consist of four isome- 
trics and the dilation algebra. In the limit of small mass of point gravitating 
source, e = 2M [12], all lost conservation laws would be recovered as first or- 
der trivial approximate symmetries. For the orbital equation, the first-order 
approximate symmetries are given in Appendix B. 

3 Bardeen Model 

Reissner-Nordstrom solution is a static, spherically symmetric and asymptot- 
ically fiat solution of the Einstein-Maxwell field equations having a singularity 
structure. There arise some problems when one applies the laws of physics at 
the singularity point (origin of the radial coordinate) of the RN solution. To 
resolve this problem, Bardeen (1968) constructed a singularity free solution 
[301 131] using the energy-momentum tensor of non-linear electrodynamics as 
the source of the field equations. The Bardeen model is known as a regu- 
lar black hole solution because for a particular ratio of mass to charge, this 
represents a black hole and a singularity free structure. After that many 
regular black hole solutions have been discussed by various authors, for ex- 
ample, [32] . These metrics, being analogues to the RN metric with magnetic 
charge, have a quite similar global but regular structure as that for the RN 
spacetime. 

3.1 Isometry Algebra and Orbital Equation of Motion 

The line element representing Bardeen model is given by [31] 

The asymptotic behavior of the solution can be observed by the simple ex- 
pansion of the term gu as follows 

2m 3me^ ^ / 1 n 
9tt = l + + O - 
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where m is the mass of the configuration. When r — > oo, the parameter 
e does not vanish as the Columbian term ^ but vanishes as the term ^. 
Therefore, one cannot recognize the parameter e as the Coulomb charge Q 
of the RN spacetime for which gu is given by 

2m g2 
9tt = l h 

Indeed, in the Bardeen model, the parameter e represents the magnetic 
charge of the non-linear self-gravitating monopole [33] while we have taken 
the gravitational units, i.e., G = c = 1. For e = 0, it turns out to be the 
Schwarzschild spacetime while for m = = e, it becomes Minkowski space- 
time. The simultaneous solution of the Killing equations for the Bardeen 
model is given by 

= Co, = 0, k"^ = c\ sin — C2 cos 0, )€" = cot 9{ci cos + C2 sin 0) + C3. 

This solution corresponds to the generators Yq, Yi, Y2 and Y3 which form 
an algebra so(3)©i? implying the conservation laws for energy and angular 
momentum only while the linear and spin angular momentum conservation 
laws are lost. We would check whether one can recover these lost conservation 
laws in the limit of small gravitational mass and charge or not. 

The set of geodesic equations for the spacetime ([9]) can be written as 

t + 2[ ^ , i P— -J-]tr = 0, 10 

1- (^,-2 _|_ g2^5/2 -1^ ~ _|_ g2^3/2-l^ 

mrir'^ — 2e^) , 2mr^ , , 9 



_|_ g2^5/2 L _|_ g2^3/2J 

--[i-pf^](^' + --'^<^') = 0' (11) 

-r^-sin^cos^02 = 0, (12) 
r 

0+ -f0 + 2cot^^0 = 0. (13) 
r 

Using the spacetime, four geodesic equations and the symmetries of the equa- 
torial plane, we obtain the following relativistic equation of motion (orbital 
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equation of motion) 

(fu 2mu'^ m(l — 2e^M^) mu'^(l — 2e^v?') 

^ ^ " " (1 + M2e2)3/2 + /^2(l + 62^2)5/2 + (1 + e%2)5/2 

where h is the classical angular momentum per unit mass and u = ^ 

3.2 Perturbation Parameters and Approximate Sym- 
metries of the Orbital and Geodesic Equations 

The Bardeen model is a singularity free solution and represents black hole 
for the following inequality 

27e^ < 16ml (15) 

Since the vanishing of mass and charge provide the Minkowski spacetime, so 
the perturbation parameters are to be defined in terms of mass and charge 
of the black hole to recover the lost conservation laws. The required pertur- 
bation parameters are 

4 

e = 2m, = ke^] < k < —, 

(where we have used the inequality defined in Eq. ffTSj) ). Introducing these 
parameters in the Bardeen model, we obtain the following third-order per- 
turbed spacetime 

= [1 - ; + 1^1*' - [1 + ; + ^ - 1^1*^ - + -"^ Oil'')- 

Retaining only the terms of first and second order in the above equation and 
neglecting O(e^), it reduces to the second-order perturbed spacetime for the 
Schwarzschild solution. Also, it reduces to the Minkowski spacetime in the 
limit of e = 0. The third-order perturbed orbital equation of motion can be 
written as 

d^u , 1 3 n, n.l'^ku^ dku^ . 

For this equation of motion, the exact and first-order parts of the sym- 
metry generator are the same as for the Schwarzschild spacetime given by 
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Eqs. ( lAlOj) - (1A16P . Since there is no term quadratic in e, therefore its second- 
order approximate symmetries are zero. Now we calculate its third-order 
order approximate symmetries. For this purpose, we apply the operator 

to Eq. fll6l) . then substitute the values of the prolongation coefficients and 
retain only the terms involving e'^, we obtain a system of four DEs. The 
simultaneous solution of these equations is exactly the same as given by 
Eq. (IA10j) - (IA14l) . Thus we can conclude that there exist no non-trivial third- 
order approximate symmetries. 

In order to recover full set of lost conservation laws, we apply the ap- 
proximate symmetry procedure to the system of geodesic equations given by 
Eqs.f lT0|) -f lT3|) . The third-order perturbed geodesic equations are 

t+-^if + ^if + ^(2-9A:)ir = 0, (17) 

f - r(9^ + sin^ e<p^) + e[^(P + r^) + 6^ 

+ sin2^</>2] - —It' - r'] - ^[— P - (1 - ^)f2 
^ ^ 2r^ ^ ^ 2r4 ^ 2 ^ 2 ' 

-3A;r2(^2^sin2^02)] =0, (18) 

.. 2 • 

e + -r^-sin^cos^(/)2 = 0, (19) 
r 

+ -f(/) + 2cot^^0 = 0. (20) 
r 

For e = 0, this system of equations reduces to the system for the Minkowski 
spacetime. For = with e 7^ 0, this system results in the first-order per- 
turbed system of geodesic equations for the Schwarzschild spacetime. There- 
fore, its exact and first-order approximate parts of the symmetry generator 
are the same as for the Minkowski and Schwarzschild spacetimes respectively. 

For the second-order approximate part of the symmetry generator, we 
apply the second-order prolongation of the symmetry generator X^^' . This is 
given by 
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where ^ = + + and 7]" = r]^ + tr]l + e^i]]^, {u = 0, 1, 2, 3). The 
apphcation of this operator to the second-order perturbed system of geodesic 
equations (i.e., = 0, e 7^ 0, 7^ 0) yields the system of 60 DEs. The 
solution (where we have retained only the terms involving e^) to this system 
of DEs is 10 Killing vectors and the dilation algebra ^2. 

For the third order approximate symmetries, we apply the above operator 
X'^] to third-order perturbed geodesic equation given by Eqs. ffT7|l - fl2n]) . where 
^ = 'Co + e6 + e^6 + e^6 and rj" = -q^ + er]'^ + e^r]^ + e^rji^ . The application of the 
operator again yields a system of 60 determining equations (where we have 
retained only the terms involving e^). In the construction of this system 
of DEs, we have used 10 Killing vectors in which 4 symmetries (given by 
Eqs. (IAip - (IA3p ) are taken as exact and the remaining 6 (given by Eqs. (IA4p )- 
(IA9p ) as first-order and second-order parts of the symmetry generators. 

We have seen that out of four constants corresponding to exact symmetry 
generators, two have been appeared and canceled out while the six constants 
corresponding to first and second-order parts of the symmetry generators 
appear. All these six constants turn out to be zero when we solve this 
system of equations simultaneously and it results in the system of DEs for 
the Minkowski spacetime. Its solution provides 10 Killing vectors and the 
generators corresponding to dilation algebra d2- This leads to no non-trivial 
symmetry in third-order approximation. For this model, the constants corre- 
sponding to dilation algebra ^0(5) = CqS + Ci do not cancel automatically but 
collect some expression for cancelation. This expression is given as follows 

1 = 1 -, (22) 

where we have used the value of k. Since ^(s) is the coefficient of ^ (where 
s corresponds to proper time) in the symmetry generator given in Eq.([2]), 
so we call this factor as energy re-scaling factor of a test particle due to the 
presence of charge and mass. 



4 Summary and Discussion 

This paper is devoted to study energy of the Bardeen model (regular black 
hole solution) by using the approximate symmetries. These solutions have 
the isometry algebra so(3)©i? while the system of geodesic equations have 
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algebra so{3)Q)R®d2- The strength of the approximate symmetries (i.e., per- 
turbation parameters) provide the information about the amount of energy 
contained in the spacetime field. The second-order approximation to symme- 
try generator yields only the lost conservation laws and hence no information 
regarding the energy re-scaling factor is obtained. For the energy re-scaling 
factor, we have to evaluate the third-order approximate symmetries. 

Firstly, we have evaluated the third-order approximate symmetries of the 
orbital equation of motion. The exact and the first-order parts of the sym- 
metry generators are the same as for the Schwarzschild spacetime. However, 
the second-order approximate part of the symmetry generator is zero. In the 
third-order approximation, one can obtain the trivial symmetries only (i.e., 
of Minkowski spacetime) given by Eqs. (lA10p - (IA14p . 

Secondly, we have calculated the third-order approximate symmetries of 
the perturbed geodesic equations. Here the exact, first-order and second- 
order approximate parts of the symmetry generator turn out to be the same 
as for the Schwarzschild spacetime. The third-order approximate symmetries 
of the third-order perturbed geodesic equations provide no non-trivial sym- 
metry generator. However, we have obtained energy re-scaling factor given 
by Eq. (!22|) . The re-scaling factor is of great interest as it shows the reduction 
in energy by the ratio of magnetic self-energy of the source to its gravita- 
tional self-energy. For e = 0, it vanishes and coincides with the Schwarzschild 
spacetime. 

When a particle is placed in the field produced by charge gravitating 
source, some force is exerted on it which should be position dependent. How- 
ever it is not clear whether the modification in mass (energy) in the presence 
of charge should be position dependent. The expression for re-scaling of force 
calculated through the pseudo-Newtonian formalism [31] is given by 



It is position dependent and provides the reduction in force by a factor de- 
pending on the ratio of the magnetic potential energy at a distance r to 
the rest energy of the gravitational source. However, it is not a convenient 
expression for relating the magnetic self-energy of the configuration to its 
gravitational self-energy. It should be position independent where the radial 
parameter r would be canceled out. Our obtained re-scaling factor is clearly 
position independent and corresponds to change in mass (energy) due to the 
presence of magnetic charge. Therefore, it provides a physically more sig- 
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nificant expression for relating these self-energies as compared to the results 
found through different approaches [33]. 

We would like to mention here that Bardeen model with magnetic charge 
is the analogue of the RN spacetime with electric charge for which, in the 
second-order approximation an energy re-scaling factor 

1 - 2A; = 1 - -i— (23) 

is obtained. However, in contrast to this, we are forced to calculate its third- 
order approximate symmetries because the second-order approximation pro- 
vide no information about the energy re-scaling factor. Clearly both these 
scaling factors are position independent providing more significant results. 
For the orbital equation of motion, the second-order approximate symme- 
tries are zero (while for the RN spacetime, there exists non-zero and trivial 
second-order symmetries). In literature p3| [16] . a difference between the 
conservation laws for the full system of geodesic equations and single orbital 
equation of motion is noted. This difference also holds for the Bardeen so- 
lution. It would be worthwhile to investigate the approximate symmetry 
generators and the concept of energy re-scaling for the rotating stringy black 
hole solutions which contain electric, magnetic and both charges. 



Appendix A 

The symmetry generators of the Minkowski spacetime are given as follows: 
Yo = ^, Yi = cos(/)^-cot0sin0^, (Al) 
Y2 = sin0— + cot^cos0— , (A2) 

OU 0(j) 

. ^ ,9 COS 9 cos d) d CSC Osmd d , , , , 

Y4 = sm^cos0— + ^— ^ — , A4 

or r 06 r 0(p 

. ^ . , (9 cos 6^ sine/) d esc 6* cos d , . ^, 

Y5 = sm0sm0— + ^— + ^ — , (A5) 

or r ou r od) 
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„ I d . ^ , d cos 6*008 d 080 6*81110 d , 
Y, = rsm^cos0- + t(sm^cos0- + ^^- ^-), 

(A7) 

• /I • , ^ / ■ n ■ , 9 cos ^ sin 9 

Yg = r8in6'sm0— + tfsmt^sm^— — I — 

ot or r o9 

CSC 7? cos d . , . 

, ^ d sm.9 d . , . , 

Yg = rcos^— + i cos^^ — . A9 

at or r aO 

Appendix B 

For the orbital equation, the first-order approximate symmetries of the Schwarzschild 
spacetime include 

d d d 

Yo = u—, Yi = cos0— , Y2 = sin(/.— , (AlO) 
ou ou ou 

Ys = ^, Y4 = cos20— -Msin2</.— , (All) 

Ys = sin 20— + u cos 20— , (A12) 
O0 ou 

d d 

Yg = MCOS0— — u^sin0— , (A13) 
ocp ou 

d d 

Yy = usin0— + M^cos0— , (A14) 
a0 ou 

and two non-trivial stable approximate symmetry generators are 

d d d 

Yai = sin0— + e(2sin0— + MCOS0— ), (A15) 

d d d 

Ya2 = COS0- e(2cos0— - wsin0— ). (A16) 

ou 0(f) ou 
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